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Abstract — The algebraic formulation for linear network coding 
in acyclic networks with the links having integer delay is well 
known. Based on this formulation, for a given set of connections 
over an arbitrary acyclic network with integer delay assumed for 
the links, the output symbols at the sink nodes, at any given time 
instant, is a -linear combination of the input symbols across 
different generations, where ¥q denotes the field over which the 
network operates. We use finite-field discrete fourier transform 
(DFT) to convert the output symbols at the sink nodes, at any 
given time instant, into a -linear combination of the input 
symbols generated during the same generation. We call this as 
transforming the acyclic network with delay into n-instantaneous 
networks {n is sufficiently large). We show that under certain 
conditions, there exists a network code satisfying sink demands 
in the usual (non-transform) approach if and only if there exists a 
network code satisfying sink demands in the transform approach. 
Furthermore, we show that the transform method (along with the 
use of alignment strategies) can be employed to achieve half the 
rate corresponding to the individual source-destination min-cut 
(which are assumed to be equal to 1) for some classes of three- 
source three-destination unicast network with delays, when the 
zero-interference conditions are not satisfied. 

I. Introduction 

Network coding was introduced in [ 1 ] as a means to improve 
the rate of transmission in networks. Linear network coding 
was introduced in [2] and it was found to be sufficient 
to achieve the maxflow-mincut capacity in certain scenarios 
such as multicast. The existence problem of network coding 
for networks without delay was converted into an algebraic 
problem in [3]. The case of acyclic networks with delays was 
abstracted in [3] as acyclic networks where each link in the 
network has an integer delay associated with it. 

The problem of network coding for multiple unicast sessions 
was considered in [4], [5]. In [6], the concept of interference 
alignment from interference channels [7] was extended to 
instantaneous unicast networks with three source-destination 
pairs for the case where, each source-destination pair has a 
min-cut of 1. This was called network alignment and it is 
useful in guaranteeing a mininum throughput when the zero- 
interference conditions in Theorem 6 of [3] cannot be satisfied. 

The motivation behind this work is striving to provide a 
minimum throughput guarantee when the zero-interference 
conditions cannot be satisfied in an acyclic network with delay, 
while not making use of any memory at the intermediate nodes 
(i.e., nodes other than the sources and sinks). The set of all Fg- 
symbols generated by the sources at any particular time instant 
are said to constitute the same generation. The output symbols 
at the sink nodes, at any given time instant, is a Fq-linear 



combination of the input symbols across different generations, 
where denotes the field over which the network operates. 
We convert the output symbols at the sink nodes, at any 
given time instant, into a F^ -linear combination of the input 
symbols generated during the same generation, by using tech- 
niques similar to Multiple Input Multiple Output-Orthogonal 
Frequency Division Multiplexing (MIMO-OFDM) [8]. We 
call this technique as the transform technique, since we use 
DFTs over finite fields towards achieving this instantaneous 
behaviour in the network. As a first step towards guaranteeing 
a minimum throughput when the zero-interference conditions 
cannot be satisfied in an acyclic network with delay, we 
consider a three-source three-destination unicast network with 
the source-destination pair denoted as Si-Di {i £ {1,2,3}). 
We also assume a min-cut of one between source Si and 
destination Di. Under this set-up, we apply the transform 
techniques and network-alignment to find conditions under 
which the network can guarantee a throughput close to half 
for every source-destination pair Si-Di {i G {1,2,3}). This 
method does not make use of memory at the intermediate 
nodes. 

The contributions of this paper are as follows. 

• We convert the output symbols at the sink nodes, at any 
given time instant, into a F^ -linear combination of the 
input symbols generated during the same generation using 
finite-field Discrete Fourier Transform (DFT). We call 
this as transforming the acyclic network with delay into 
n-instantaneous networks, where, n is sufficiently large. 

• Using a constructive proof, we show that there exists a 
network code (satisfying a certain property) that achieves 
the sink demands in the usual (non-transform) approach 
if and only if there exists a network code satisfying sink 
demands in the transform approach . 

• For a three source-three destination unicast network with 
delays, which do not satisfy the zero-interference con- 
ditions, we extend the transform techniques to achieve 
atleast half the rate corresponding to the individual 
source-destination min-cut (which are assumed to be 
equal to 1), along with the use of alignment strategies. 
In particular, the contributions for the three source-three 
destination unicast network with delays are as follows. 

1) When the min-cut between Si-Dj is greater than 
or equal one, V € {1,2,3} (z ^ j), we 

derive sufficient conditions under which network 
alignment can achieve half the rate corresponding 



to the individual source-destination min-cut, with 
time-invariant Local Encoding Kernels (LEKs). 

2) The network alignment procedure with time- 
invariant LEKs is then generalized with the use of 
time-varying LEKs. 

3) When the min-cut between Si-Dj is zero for some 

G {1,2,3} {i ^ j), we derive sufficient con- 
ditions under which network alignment can achieve 
atleast half the rate corresponding to the individual 
source-destination min-cut. 

The organization of this paper is as follows. In Section II, 
we review the system model for acyclic networks with delays 
presented in [3]. Section III presents the central contribution 
of this work, i.e., the transform technique using which we 
convert the usual convolutional behaviour of the network 
into instantaneous behaviour. In Section HI, we also prove 
the interchangeability of solving the usual (non-transform) 
network code existence problem and the counterpart in the 
transform technique. In Section IV, we combine our trans- 
form technique with the alignment techniques for acyclic 
instantaneous networks given in [6] to achieve an asymptotic 
throughput of 1/2 for certain classes of acyclic networks with 
delays, even when the zero-interference conditions cannot be 
satisfied in such networks. We conclude our paper in Section 
V with a discussion and directions for further research. 

Notations: The cardinality of a set E is denoted by \E\. 
A superscript of t accompanying any variable (for example, 
e^*)) or any matrix (for example, M^*)) denotes that they are a 
function of time t. The row, column element of a matrix 
A is denoted by [A]ij. The notation P C Q denotes that the 
columns of the matrix P are a subset of the columns of the 
matrix Q. Span(_P) indicates the sub-space spanned by the 
columns of the matrix P. The determinant of a square matrix 
A is denoted by det{A). An identity matrix of size /i x /i 
is denoted by 7^. For three-source three-destination unicast 
networks we shall use the term destination to denote sink. A 
Galois Field of cardinality is denoted by GF{p™) where, 
p is a prime number and m is a positive integer. 

II. System Model 

First, we shall briefly review the system model from [3]. We 
consider a network represented by a Directed Acyclic Graph 
(DAG) g = {V,E), where V is the set of nodes and E is 
the set of directed links. We assume that every directed hnk 
between a pair of nodes represents an error-free Unk and has a 
capacity of one F^, symbol per link-use. Multiple links between 
two nodes are allowed and the i*'' directed link from vi G V 
to ^2 G 1^ is denoted by {vi,V2, i). The head and tail of a link 
e = (t'l. ('2, i) are denoted by V2 = head(e) and vi = tail(e). 
A link between a pair of nodes can have an arbitrary finite 
integer delay. Let X{v) = {X{v, l),X{v, 2), X{v, be 
the collection of discrete random processes that are generated 
at the node v. Let Xy_ = [X{v,l) X{v,2) ... X(i>,/z„)]^. 
The random process transmitted through Unk e is denoted by 
Z{e). Communication is to be established between selected 
nodes in the network, i.e., we are required to replicate a subset 



of the random process in X{v) at some different node v'. A 
connection c is defined as a triple {v,v' ,X{v,v')) € V x 
V X Vx(v)^ where Vx{v) denotes the power-set of X{v). For 
the connection c, v is called the source and v' is called the 
sink of c, i.e., v = source(c) and v' = sink(c) (source(c) 
7^ sink(c)). The collection of Vy' random processes y{v') = 
{Y{v' ,l),Y{v' ,2), ...,Y{v' ,1/^,)} denotes the output at sink 
v'. Let = [Y{v', 1) Y{v', 2) ... Y{v', 

The input random processes X{v,i), output random pro- 
cesses Y{u,j) and random processes Z{e) transmitted on 
the link e are considered as a power series in a delay 
parameter D, i.e., Xiv,i) = YT=o^''*\'",i)D\ Y{u,j) = 
J:ZoY^'Hu,j)D\ and Z{e) = Z^^) {e)DK 

Let 5 = (y, £■) be an acychc network with arbitrary finite 
integer delay on its links. ^ is a F^-Unear network [3], if for 
all links the random process Z{e) on a link e = {v, u,i) £ E 
satisfies 

Z(*+i)(e) = £a,,eXW(^;,i) + 0e',eZ^'\e') 

j—l e':head(e')— tail(e) 

where, aj^e and (3e',e belong to F,, where q = p™, for some 
prime number p and positive integer m > O.The output at any 
sink node v', is taken to be 

Y('+'\v',j)= Yl ^e'jZ('\e') (1) 

e':head(e')=i'' 

where eg'.j £ F,. The coefficients, otj^e, Pe',e and e^'j are also 
called local encoding kernels (LEKs). The vector consisting of 
all LEKs is denoted by e. Note that in [3], the definition for the 
output processes at any given time instant at any sink involves 
linear combinations of the received processes and output 
processes across different previous time instants, and hence 
the variables involved in such linear combinations together 
performed the function of decoding the received processes at 
the sinks to the demanded input processes. However, in (1), 
at every sink, we only define a preprocessing of the received 
symbols corresponding to the previous time instant alone. The 
outputs as t varies, will then be used by sink-j 

to decode the demanded input processes using sufficient delay 
elements for feed-forward and feedback operations. These 
LEKs are time-invariant unless mentioned otherwise. 

We assume some ordering among the sources so that the 
random process generated by the sources can be denoted, 
without loss of generality, as XjJ^D), X2{D), XJ^D), 
where ,s denotes the number of sources and XjiD) is a /i^ x 1 
column vector given by 

)U{D) = [X^l{D) Xa{D) ...X^^^{D)f. 

Similarly, we assume some ordering among the sinks so that 
the output random process at the sinks can be denoted, without 

loss of generality, as Yi{D), Y2{D), Yr{D), where r 
denotes the number of sinks and YJ^D) is a t-j x 1 column 
vector given by 

Yd,D) = [Y,i{D)Y,2{D) ...Y,,^{D)f. 



Let 



Y{D) = [Y,{Df Y2{Df ... YriDff 
= [yi{D) y2iD) ... y.{D)f . 



Also, let 



X{D) = [X^{Df X2{Df ... Xs^^Dff 
= [.Tip) X2iD) ... x^iD)f, 

where /x = J2i=i f^i ^ = Si=i Henceforth, the tail 
of an edge originating from a source will be identified by the 
source number and the head of an edge terminating at a sink 
will be identified by the sink number. From [3], we have 



Y{D) = M{D)X(D) 



(2) 



where, M{D) denotes the network transfer matrix of size 
V X jjL with elements from Fg[£'], the ring of polynomials in 
variable D with coefficients from F,. Now, M{D) can also 



be written as 



M{D) 



Mn(D) 
Mi2{D) 



M2i{D) 

M22{D) 



Mir{D) M2r{D) 



M,i{D) 
M,2{D) 



(3) 



where My {D) denote the network transfer matrix from source 
i to sink j and is of size Vj x /Xj. Let d^^j^ and denote the 
maximum and the minimum of all path delays from source-i 
to sink-j, V between which a path exists. Let 



dmax — '^macc 

Hence, M{D) can be written as 



M{D) = 



E 



. d=0 / 



where M^''^ G F^^^ represents the matrix-coefficients of D"^ 
of the polynomial elements of M(D). 

Since Z)'^™'" just adds a constant additional delay to all the 
outputs, we can take, without loss of generality, M{D) as 



Hence, (D) can be alternatively written as 

Mij{D) 



d=0 



(4) 



(5) 



For each sink-j, we also define Mj{D) to be the //, x /i 
submatrix of M[D) that captures the transfer function between 
all the sources and the sink-j, i.e., 

M,{D) = [M^,{D) M2,{D) ... M,,{D)]. (6) 

In the network Q, let Cj denote the set of all connections to 
sink-j. Let C — U^^^Cj. The following lemma from [3] states 
the conditions for solvability of acylic networks with delay. 



Lemma 1 ( [3]): An acycUc network with delay is solvable 
iff there exists an assignment to the LEKs e such that the 
following conditions are satisfied. 

1) Zero-Interference: Ml^\li) = 0, for all pairs (source-i, 
sink-j) of nodes such that (source-i, sink-j, Xi^^^\D)) 
^ Cj for all < d < dmax, where Ml'^\li) denotes 
the If column of M^f and Xi^^'\D) denotes the if 
element of Xi{D). 

2) Invertibility: For every sink-j, the square submatrix 
M'j{D) of Mj{D) formed by juxtaposition of the 
columns of Mij{D) (V i; I < i < s) other than those 
involved in the zero-interference conditions is invertible 
over ¥q{D), the field of rationals over Fg. 

A network code for {G, C) is defined to be & feasible network 
code if it achieves the given set of demands at the sinks 
i.e., if the zero-interference and the invertibility conditions are 
satisfied. 

A. System Model for time-varying LEKs 

When the LEKs are time-varying, we can't express the 
input-output relation as in (2). Hence, first, we need to derive 
the input-output relation involving transfer matrices which are 
dependent on varying LEKs. Retaining the notations as already 
introduced, we only point out the changes in the system model 
here. 

Without loss of generality, we assume that a link between a 
pair of nodes has a unit delay (if the link has any other non- 
zero integer delay, we could introduce an appropriate number 
of dummy nodes in between the pair of nodes which are then 
connected by links of unit delays). For a given DAG Q with 
integer delay on its links, define the adjacency matrix of Q at 
time t as the x |£| matrix whose elements are given 
by 

= / ^Se,- head(ei) = tail(ej) 
^ \ otherwise 

Let the entries of ijl'><\E\ matrix at time t, be given by 



1 (t) 



X 

otherwise 



Also, let the entries of vx \E\ matrix B^*\ at time t, be given 
by 

r it) ^ 

otherwise 



[B 



(Ol.. - 



Let the set of vectors denoted by ej^i-t^) j-jjg denote the set 
of LEKs from time instant ti to time instant t2 (t2 > ii), i.e., 

where denotes the LEKs at time ti. Since the LEKs are 
time varying, the network transfer matrix is given by 

M(D,t)T'=(A(*-i)7i>fA(«-2);<r(*-i)D2+A(*-3);<r(*-2)K(t-i)D3 



where M^^^ = 0, i.e., the zero matrix, as each link in the 
network is assumed to have a unit delay. 

Since acyclic networks with delay are analogous to multiple 
transmitter-multiple receiver Unear channel with time-varying 
impulse response between every transmitter and every receiver, 
the output symbols for the acychc network with delay, at time 
instant t, at sink-j, is given by 



-d) 



(7) 



i=l d=0 



III. Transform Techniques for Acyclic Networks 

WITH Delay 

In this section, we show that the output symbols at all the 
sinks which was originally a Fg-linear combination of the 
input symbols across the different generations, at any given 
time instant, can be transformed into a Fg-hnear combination 
of the input symbols across the same generation. 

Consider a matrix A of size nv x nn given by 



■^0 




Al-1 

Al-2 

Al 



Al 
Al-1 







Al 







where, AiS (0 < z < L) are matrices of size u x n, whose 
elements belong to Fg and n >> L. Note that the [i + 1)* 
row of matrices is a circular shift of the row of matrices in 
A. We assume that n divides g — 1. Since q = p™, p and n are 
coprime. The choice of n is such that, there exists an a e Fg 
such that n is the smallest integer for which a" = 1. This is 
indeed possible [9]. Define matrices Aj {0 < j < n — 1), of 
size u X fi, as 



i=0 



Let F be the finite-field DFT matrix given by 



1 1 

1 Q 

1 a 



Also, define the matrix as 















alij, 














.1^ 







a 

^2(n-l) 



^(n-l)(n-l) 



I, 



(8) 



Similarly we can define matrix Q^. The following theorem 
will be useful in establishing the results subsequently. 
Theorem 1: The matrix A can be block diagonalized as 

A = QiiAQ^ , 



where, A is given by 

in-i 

i„-2 







... ■ 

... 

io. 



Proof: Proof is given in Appendix A. ■ 
Now, consider an arbitrary acychc network with delay. From 
(2) and (3), 



^iD) = J2Mij{D)Xi{D). 



(9) 



Now, consider a transmission scheme, where we take n 
(>> dmax) generations of input symbols at each source and 
first transmit last d^ax generations (which we call the cyclic 
prefix) followed by the n generations of input symbols. Hence, 
n + dmax time slots at each source are used for transmitting 
n generations. Then, (9) can be written as (10) using (4). 
Now, after discarding first dmax outputs at sink j, (10) can be 
re- written as (11) (given at the top of the next page). Using 
Theorem 1, (11) can be re- written as 



i=l 



(12) 



where. 



y.in-ir 

v.("- 

^3 



2) 



Mi. 



M, 



y.(0) 











p^(n-l)n 






(0) 



M. 



Qn^Xi"- instead 



Now, at each source i, transmit X^'" 
of Xi"'. Then, at each sink j, we receive 1^'". Let 1^" = 
Q-ir/". Then, from (12), ~ ~ 

s 

s 

Yi'=Y,MijXr (13) 



i=l 



Now, (13) can be re-written as (for <t <n — 1) 



(14) 



Hence, each element of Y,^*^ is a F^ -linear combination of 
the input symbols across the same generation. We now say 



FM. 



E 



(0) 



(1) 



M, 



(0) 



(1) 



M, 



(0) 



Mi: 



My 

(d, 
m;, ' 



MK 



(10) 













= E 

■i=i 







M}, 



(1) 



M^: 



M}: 



-1) 

-2) 



M, 








M}" 



(0) 



(11) 



that we have transformed the aycUc network with delay into 
n-instantaneous networks. 

Remark 1: Note that the linear processing of multiplying 
by matrices Q^^ at source-i and at sink-y are done in 
a distributed fashion which is necessary because the sources 
and sinks are distributed in the actual network. 

Remark 2: One can observe that transmitting X'^^ = 
QixiXi" implies taking DFT across n generations of each 
of the Hi random-processes generated at source- i. Similarly, 
the pre-multiplication by Ql^^ at sink-j simply implies taking 
IDFT across n generations of each of the vj random-processes 
received. The entire processing, including addition of cyclic 
prefix at source-i and removal of cyclic prefix at sink-j is 
shown in a block diagram in Fig. 1 (given at the top of the 
next page). 

Now, let us re-write (14) as 

i=l li = l 

where Ml-\li) denotes the if column of M^*^ and Xi^^^k) 
denotes the if element of 

Similar to the zero-interference and invertibility conditions 
in Lemma 1, we have the following theorem for solvabiUty of 
(14). 

Theorem 2: An acylic network {Q, C) with delay, incorpo- 
rating the transform techniques, is solvable iff there exists an 
assignment to variables e such that: 

1) Zero-Interference: Ml^\li) = for all pairs (source-i, 
sink-j) of nodes such that (source-i, sink-j, X^^\li)) 
^ for < t < n - 1. 

2) Invertibility: If Cj contains the connections (source-ii, 
sink-j, Xi^*\k^)), (source-i2, sink-j, Xi^''\U^)), ■ • ■ , 
(source-iss sink-j, Xi^,^^\li^,)), then, the sub-matrix 

M^l^jiK,)] is a nonsingular 
I'j X i^j matrix for < t < n — 1. 



The network code which satisfies the invertibility and the 
zero-interference conditions for {G,C) in the transform ap- 
proach using a suitable choice of a for the DFT operations is 
defined as an feasible transform network code for {Q,C). 

A. Existence of a network code in the transform approach 

In this section, we prove that under certain conditions there 
exists a feasible network code for a given {Q, C) if and only if 
there exists a feasible transform network code. Towards that 
end, we prove a lemma. We first define the polynonaial f{D) 
which will be used henceforth throughout this paper. 

r 

f{D) = Y{det{M'^{D)). (15) 

j=i 

where, Mj{D) is the square submatrix of Mj(D) indicating 
the source processes that are demanded by sink-j'. 

Lemma 2: Suppose there exists a feasible network code for 
iG,C) over some field F^. For some a € ¥qa (for some 
positive integer a), the local encoding kernels defined by the 
feasible network code for {Q, C) (viewed in the extension field 
Fgo) along with the DFT operations defined using a result in 
a feasible transform network code for {Q,C) if and only if 
/(a*) ^ for all < i < n - 1. 

Proof: Proof is given in Appendix C. ■ 

We now prove the following theorem which concerns with 
the relationship between the existence of a feasible network 
code and a feasible transform network code for {G,C)- 

Theorem 3: Let {G,C) be the given acyclic delay network 
with the set of connections C demanded by the sinks. There 
exists a feasible transform network code for {Q, C) if and only 
if there exists a feasible network code for {G,C) such that 
{D-l)\f{D). 

Proof: Proof is given in Appendix D. ■ 

Remark 3: Based on the constructive proof of Theorem 3, 
a large field might be required for the existence of a suitable 
value for a that defines the necessary transform for the net- 
work, under the condition that the rate-loss ( ) due to the 
transform approach be less. The transformed network would 
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Fig. \. Block Diagram to illustrate linear processing at Suurcc-i and Sink-;. 
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then have to be operated over this large field, i.e., the matrices 

' (t) 

Mij have elements from this large field (which is at least 
a degree n extension over the base field over which the non- 
transform network code is defined). It is known that (see [11], 
for example) inverting a Vj x matrix (at some sink-j) takes 
0(z/?) computations, however over the extension field. In the 
process of computing these inverses, the information symbols 
corresponding to the n generations are obtained by Gauss- 
Jordan elimination. In terms of base field computational com- 
plexity, the complexity of computing the inverse of the transfer 
matrix becomes O (^'?n(logn)(loglogn)), as each multipli- 
cation in the extension field involves O (n(logn)(loglogn)) 
computations over the base field [12] (it is equivalent to mul- 
tiplying two polynomials of degree at least n — 1 over the base 
field). The total complexity of recovering the input symbols 
at all the n generations is then O (n^t'|(logn)(loglogn)). 

On the other hand, if the non-transform network code 
is used as such, the transfer matrices Mj(D) consist of 
polynomials of degree upto d„,„ r in D over the base field. 
Again, it is known (see [11], for example) that finding the 
inverse of such a matrix has complexity 0{iyjdma:i)- To do a 
fair comparison with the transform case, we consider decoding 
of n-generations (n being large as in the transform case) of 
information. Note that inversion of the matrix M'j{D) does 
not give us the information polynomials directly. A naive 
method of obtaining the each information polynomial would 
then require multiplications of polynomials over the base 
field (each of which has complexity O (n(logn)(loglogn)), 
assuming that Vjdmax < n.) and one division between polyno- 



mials (again with complexity O (n(logn)(loglogn))). There- 
fore, the total complexity involved in recovering the infor- 
mation sequences would then be O (t'^n(logn)(loglogn)) -|- 
O (i/jn(logn)(loglogn)) + 0{i'^dmax) computations. 

Thus, we see that there is an advantage in the complexity 
of decoding in the non-transform network compared to the 
transform network (inspite of using the least possible size for 
the extension field). Therefore, complexity reduction is not an 
advantage of the transform process. 

We now present an example acyclic network in which there 
exists a feasible network code, using which we obtain a 
feasible transform network code for some choice of n > 7. 

Example 1: Consider the network Q shown in Fig. 2. This 
is a unit-delay network (where each edges have a delay of one 
unit associated with it) taken from [13]. For 1 < i < 3, each 
source .s^ has an information sequence Xi{D). This network 



has non-multicast demands, with sinks Uj 



1 < 



< 3 



requiring all three information sequences, while sink m re- 
quires {xi{D) ^ xj,{D)} and demands {x2{D),xt,{D)} . 
Let C denote these set of demands. A feasible network code 
for {Q,C) over F2 as obtained in [13] can be obtained by 
using 1 as the local encoding kernel coefficient at all non-sink 
nodes. The transfer matrix M„^(£)), the invertible submatrix 
M^.{D) of Muj{D), and their determinants for the sinks 
Uj : 1 < .7 < 5 are tabulated in Table 1 . 

We therefore have /(£>) = £)25 ^^^^ (^at /(I) ^ and 
dmax = 4 for this network. Therefore, with n = 2™ — 1 for 
any positive integer m > 3, i.e., a being the primitive element 
of , we will then have /(a*) ^ for any {) <t <n — 1. 



TABLE I 



Sink 


Network transfer 
matrix M^. (D) 


Invertible submatrix 
M^^^iD) ofM„^.(r>) 


Determinant of {D), 
det{M^^{D)) 


Ui 




(DO ^ 
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Fig. 2. A unit-delay network with 3 sources and 5 sinks 
By Lemma 2, we then have a feasible transfonn network code 

for {g,c). 

In the next section we shall apply these transform techniques 
to three-source three-destination unicast network with delays. 

IV. Three Source-Three Destination Unicast 
Network with Delays 

In [6], the concept of interference alignment from inter- 
ference channels [7], was extended to instantaneous acyclic 
unicast networks with three source-destination pairs for the 
case where, each source-destination pair has a min-cut of 1 
and where, the zero-interference conditions in Theorem 6 of 
[3] cannot be satisfied. This was called network alignment. 
It was shown, in [6], that for a class of such networks, it is 
possible to achieve a throughput close to 1/2 for every source- 
destination pair via network alignment. 

In this section, we deal with acyclic three source-three 



destination unicast network with delays, with each source- 
destination pair having a min-cut of 1. We employ the results 
from the previous section and show that, even when the 
zero-interference conditions of Lemma 1 cannot be satisfied, 
for a class of three source-three destination unicast networks 
with arbitrary integer delays on its links, we can achieve a 
throughput close to 1/2 for every source-destination pair by 
making use of network alignment. We take two approaches 
in achieving this - using time-invariant LEKs and using time- 
varying LEKs. 

Let the random process injected into the network by source 
Si {i e (1, 2, 3)) be Xi{D). Source Si needs to communicate 
only with destination Di (i G (1,2,3)). Here, Hi = 1 and 
i/j = l ((i,i) € 1,2,3). 

We shall consider the following two cases separately. 

1) The min-cut between source-i and sink-j is greater than 
or equal to 1, for all i ^ j. 

2) The min-cut between source-z and sink-j is equal to 0, 
for some i 7^ j. 

Case 1: The min-cut between source-i and sink-j is greater 
than or equal to 1, for all i ^ j. 

A. Achieving a Throughput of 1/2 with Time-Invariant LEKs 

Now, consider a transmission scheme, where we take 2n+l 
(>> dmax) generations of input symbols at each source and 
first transmit the last dmax generations (i.e., the cyclic prefix) 
followed by the 2n + 1 generations of input symbols. Let 
QiX?"'^^ be the input symbols transmitted by source i, where, 

j^2n+l _ rj^(2n) j^(2n-l) _ _ ^ 

Also, let ^ = ^2^2"' and 

= VgX^", where, Vi is a (2n + 1) x (n + 1) matrix. 



V2 is a (2n + 1) x n matrix, V3 is a (2n + 1) x n matrix, and 



^2 — 1^2 ^2 



■•■ X. 



/("-l)lT 
2 J 

/("-l)lT 



The quantities X^""*"^, Xg" and X3" denote the (n + 1), n, 
and n independent input symbols generated by sources- 1, 2 
and 3 respectively. Now, from (13), for j e {1, 2, 3}, 

= Myyixr+' + M2,y2xr + M3,-i^3^r> 

where, Yj^'"'^^ denotes the (2n + 1) output symbols at sink- 
j. The objective is to recover the {n + 1) independent input 
symbols of source-1, n independent input symbols of source-2 
and n independent input symbols of source-3 at sinks-1, 2 and 
3 from Y^^''+\ ^2^"+^ and Y^''+'^ respectively. 

For aycUc networks without delay, the network aUgnment 
concept in [6] involved varying LEKs at every time instant. But 
with delays it is possible, in some cases, to achieve network 
alignment even with time-inviiriant LEKs. This is what we 
show in this sub-section. 

First, note that the elements of MijS are functions of e. 

Lemma 3: Determinant of the matrix Mij V G 
{1, 2, 3} is a non-zero polynomial in e. 

Proof: Proof is given in Appendix E ■ 

Let 

= M^^M^^^^Mi^^ (k e {0, 1, .., 2n}) 
6^=) = Mi^M^fMi^ {k e {0, 1,.., 2n}) 
T = M2iM32Mi3MsiM23Mr2^ (16) 
R = MiaM^s' (17) 
S = Mi2M^2- (18) 
Now, choose 

Vi = [W TW T'^W ■ ■ ■ T"W] (19) 
V2 = [RW RTW RT^W ■ ■ ■ RT'^-^W] (20) 
V3 = [STW ST'^W ■ ■ ■ ST^'W] (21) 

where, = [1 1 • ■ • 1]"^ (all ones vector of size (2n+l) x 1). 
Since the transform approach requires that 2n + l\p"^ — 1, we 
shall find it useful in stating the exact relationship between 
2n + 1 and p which will be used in the result that follows. 

Lemma 4: The positive integer 2n -|- 1 divides p"* — 1 for 
some positive integer m iff p 1 2n -|- 1. 

Proof: Proof is given in Appendix F. ■ 

Theorem 4: The input symbols X{"^^, X!^ ,and Xg" can 
be exactly recovered at Ti, T2, and T3 from the output symbols 
y^2n+i^ if and Yi"+^ respectively subject to p f 2n -|- 1, 
if the following conditions hold. 

Rank[yi M^^M2iV2] = 2n + 1 (22) 
Rank[Mf2^M22V2 Vi] = 2n + 1 (23) 
Rank[Mi"^iM33V3 Vi] =2n+l (24) 



Proof: Proof is given in Appendix G ■ 

When the conditions of the above Theorem are satisfied, 
we say that network alignment is feasible. When network 
alignment is feasible, throughputs of ^2n+l) ' (2to+i) ' 
(on^i^ are achieved for the source-destination pairs S\ — Di, 
02 — D2, and 5*3 — D3 respectively. Hence, as n — >■ 00, a 
throughput of 1 /2 is achieved for every source-destination pair. 

Remark 4: To satisfy (22)-(24), we have to first ensure 
that Vi is full-rank. Note that 



„(2n) 









Now, any collection of [n + 1) rows of Vi is a Vander- 
monde matrix whose determinant is a non-zero polynomial 
iff VM) ^ yMT is satisfied for every fci, k2 G {0, 1, .., 2n} 
and k\^k2. But, for all the columns of V\ to be independent, 
it is enough if there exists atleast (n + 1) linearly independent 
rows. This condition is satisfied if there are atleast (n + 1) 
distinct f^s (fci G {0, 1, .., 2n}). If Vi is full-rank, then, 
by the choice of V2 and V3 in (20) and (21) respectively, 
Mf/ 71^21^2, M^2^22V2 and M^^Mz^V^, are also full-rank 
matrices. Hence, when there are atleast (n -|- 1) distinct 
{ki e {0, 1,.., 2n}), the choice of Vi, V2 and V3 atleast 
ensures that Vi, V2 and V3 are individually full-rank matrices. 

Remark 5: Note that, for three-source three-destination uni- 
cast network without delay, considered in [6], it was not 
possible to achieve network aUgnment without changing the 
LEKs with time. When there is no delay, the matrices T, R, 
and S, given in (16)-(18), would simply be equal to f{e)l2n+i 
(where, /(e) is some polynomial in e) and hence, the matrices 
Vi, V2 and as given in (19)-(21) are themselves not full- 
rank matrices. Hence, e was varied with time in [6]. However, 
in the case of delay it is easy to see from the structure of the 
matrix Mij that the matrices T, R, and S are not necessarily 
scaled identity matrices. 

The following example, taken from [6] (but with delays), 
illustrates the existence of a network where network aligrmient 
is feasible with time-invariant LEKs. 

Example 2: Consider the network shown in Fig. 3 (at the 
top of the next page). Each link is taken to have unit-delay. 
In accordance with the LEKs denoted as in the figure, the 




transfer matrices (D) are as given below. 



Mn{D) 


= M[fD^ 


= apD^, 








= M[fD^ 


+ m[1'>d^ 


= uD^ 


+ atD^ 




= mH^d^ 


= arD^ 






M2i{D) 


= M^fo^ 


= hpD\ 






M22{D) 




= UD^ 






M23iD) 


= m!§d^ 


+ a4'^d^ 


^ sD^ 


+ brD^ 


M3i(Z?) 


= M^f>D^ 






+ cpD^ 




= Mi^D' 


= ctD^, 








= Mi^D^ 


= crD^. 







Note that the network does not satisfy the zero-interference 
conditions of Lemma 1. Here, dmax = 2. It can be verified that 
network alignment is feasible with 2n + 1 = 7. Specifically, 
network alignment works with the following assigimient to the 
LEKs. 

a=b=c=p=r=t=l 

q = l + P + P^ 
n = 1 + 

where, (3 is the primitive element of GF{2^), i.e., root of the 
minimal polynomial (1 + x + x^). 

B. Achieving a Throughput of 1/2 with Time-Varying LEKs 

In this section, we shall generalize the selection of matrices 
Vt , V2 and Vz and hence Theorem 4 along with the use of time- 
varying LEKs. Throughout the sub-section we shall assume 
that the LEKs and the other variables that we shall encounter 



belong to the algebraic closure of the field Fp which is denoted 
by Fp. Clearly, once an assignment to the LEKs and variables 
are made, they belong to a finite extension of Fp. 

In this case of time-varying LEKs, the network cannot be 
decomposed into (2n + 1) instantaneous networks using the 
transform method. This is explained below. 

Consider a transmission scheme, where we take 2n + l (>> 
dmax) generations of input symbols at each source and first 
transmit last dmax generations (i.e., the cycUc prefix) followed 
by the 2n + l generations of input symbols. Let Xf^'^^ be the 
input symbol transmitted by source-?, where. 

Also, let = +1 = y2^r, and 

= Fs^r. where, Vx is a {2n+l) x{n + l) matrix, 
V2 is a (2n -\-l)xn matrix, V3 is a (2n -|- 1) x n matrix, and 

V'" — r v'C) v/("-i)iT 

^2 — 1^2 ^2 ' ' ' ^2 \ 
^3 — 1^3 ^3 ■ ■ ■ -^3 J • 

The quantities X[^'^^ , X'^ and Xg" denote the [n + 1), n, 
and n independent input symbols generated by sources- 1, 2, 
and 3 respectively. Now, from (7) and following the same steps 
involved in writing (10) and (11), for j €: {1, 2, 3}, we get (25) 
(given at the top of the next page). In brief, for j G {1,2,3}, 
we have 

= My + M2jV2X'2 + M^jV^X'^ 

where, Yj'^"~^^ denotes the {2n + 1) output symbols at sink-j 
and Mij is as given in (25). The objective is to recover the 
{n + 1) independent input symbols of source- 1, n independent 
input symbols of source- 2 and n independent input symbols of 
source-3 at sinks-1, 2 and 3 from Y^''+\ Y^"+'^ and Yi''+'^ 
respectively. 



Mff(£(2''-i'2''-i)) 



MY. 



(£!mo»-l)(^(2»-<im<,x + l,2n) 



) M~>{e 



-2)^^(2»-<im<,x + l,2n-l)-) ^(<i„ 



l)(£(2»i-<imoa!,2n-l)-j 



(25) 



(2). (-2,0)n 



M, 





jV^|ci™ax)(^(2„-l-<i„,,,2„-l)-, 



(2„) ;i,.(2n-l) 



(0)1 



Note that Mi^ is not a circulant matrix and cannot be diago- 
nalized in general. Let e' = ^£i-d.,n„a.) ^(-d^^^+i) ^ ■■,e^^"^}- 

Lemma 5: Determinant of the matrix M^j V G 
{1,2,3} is a non-zero polynomial in e'. 

Proof: Proof is given in Appendix H ■ 

Hence, the inverse of My exists. Now, let the elements of 
V\ be given by 

[Vi]ii=eij; ie{l,2,..,2n+l},jG{l,2,..,n+l}, (26) 

where % is a variable that takes values from Fp. Also, let 

V2 = M^^Mi^ViA and Vs = M^^Mi2ViB {11) 

where, the elements of the matrices A and B, of size (n + 
1) X n, are given by [A]ij = Uij and [B] 

ij ~ ^ij respectively 
{Uij and bij are variables that take values from Fp). Let 

Also, let^= {0^j\i e {l,2,..,2n + l}, j £ {1, 2, .., n + 1}}, 
a = {aij\i G {1, 2, .., n + 1}, j G {1, 2, .., n}} and & = 
{bij \i G {1, 2, .., n + 1}, j G {1, 2, .., n}}. Let 

/i {0, i', a) = det([Fi Mf/Mai^s]) 
/2(^,£',a)=dei([Mf2'M22V^2 Vi]) 
h{0,e',b)=deti[Mr^'M33V3 V^]) 
hii) = n det{Mij) 

(»j)G{l,2,3} 

f{0,e!,a,b) = h{0,e',a)f2{0,e:,a}Me,e:,b)U{e!). 

Denote the elements of a matrix C, of size n x n, by 
[C]ij = Cij, where Cij is a variable that takes values from 
% (i e {l,2,...,n},j G {l,2,...,n}). Let c = {cij\i G 
{1, 2, .., n}, j G {1, 2, .., n}}. For i G {1, 2, .., 2n + 1} and 
i G {l,2,..,n}, let 

9ij{0,e',a,b,c) = [TiV-iA]ij - [ViBC]ij. 

Let 5 ■ "''^ {9, e',a,b,c) and f;!^''''^ e' , a, 6, c) respectively de- 
note the numerator and denominator of the rational-polynomial 
gijii,£',a,b,c) {i G {l,2,..,2n+ 1} and j G {l,2,..,n}). 



Similar notation is used for the numerator and denominator 
of f {(Ij s' T k) ■ Denote f {0_, e' , a, b) and gij{d,£!,a,b,c) by 
/ and gij for short. Similar notation is used for the numerator 
and denominator of the respective rational polynomials. 

Theorem 5: For an acyclic three-source three-destination 
unicast network with delays, the input symbols X[^^'^^, 
X!^, and X3" can be exactly recovered at the sinks- 1, 
2, and 3 from the output symbols Y^''+^, and 
y^2n+i respectively, if the ideal generated by the polyno- 
mials fifj"''^ (i G {l,2,..,2n-|-l} and j G {l,2,..,n}), and 
{l - 5fi^r)fi<ir) Yl^_^ g^p^ does not include 1, where 6 is 
a variable that takes value from Fp. 

Proof: Proof is given in Appendix 1 ■ 

When the conditions of the above Theorem are satisfied, 
we say that network alignment is feasible. When network 
alignment is feasible, throughputs of (^2n+i) ' (2n+i) ' ™^ 
(2X^-1) achieved for the source-destination pairs Si — Di, 
S2 — D2, and 53 — D3 respectively. Hence, as n — >^ 00, a 
throughput of 1 /2 is achieved for every source-destination pair. 

Remark 6: If f{0,s',a,b) has to be a non-zero polynomial 
firstly, Vi has to be a full-rank matrix. This is true from the 
choice of ^1. Also, ^^"^^21^2, M-^^M22V2 and M'^M^sYs 
should also be full-rank. Since MjjS are invertible, it is 
equivalent to checking if V2 and V3 are full-rank. This is also 
true because Vi is a full-rank matrix and by choosing A and 
B as matrices that select the first n columns of Vi and the 
last n columns of Vi respectively, V2 and V3 become full-rank. 
Hence, the determinants of all the n x n sub-matrices of V2 
and V3 are non-zero polynomials. So, we have atleast ensured 
that by proper choice of Vi, V2 and V3, they are full-rank 
matrices. 

Remark 7: Note that the network aUgnment matrices in 
Section IV-A can be derived as a special case of the network 
alignment matrices in Section IV-B. Hence, Theorem 4 can 
be derived as a special case of Theorem 5. This is explained 
below. Choose eC--^--^) = gC-rfmax+i) = = ^{2n) ^ ^_ 
Also, choose the variables 9ij such that Vi in (26) takes 
the form of Vi in (19). Choose A and B, respectively, to 



be selection matrices which select the first n columns and 
last n columns of the matrices pre-multiplying them. Let 
C = In- Now, it is easy to see that TiViA — ViB is equal to 
Qi{TViA - ViB) = 0. Now, it can also be easily seen that 
the full-rank conditions in Theorem 4 are the same as saying 
that the ideal generated by (1 — 6h{9,e',a,b,c)f{9,e',a,b)) 
should not include 1. 

Case 2: The min-cut between source-i and sink-j is equal 
to 0, for some i 7^ j. 

In this case, we have totally 63 possibilities in which one 
of them is a zero-interference possibility (i.e. min-cut between 
Si-Dj is equal to zero for all i ^ j). Clearly, we need not 
consider the zero-interference possibility. 

We broadly classify the different possibiUties into four 
categories as given in Table II. All the other possibilities 
involve either permutations of the sources or require minor 
modifications in the network ahgnment procedure for one of 
these categories. For all these categories, network alignment 
can be done with time-varying LEKs too. But, the only dif- 
ference with respect to network alignment with time-invariant 
LEKs would be that the network transfer matrices cannot be 
diagonalized. 

We shall present network alignment for the categories given 
in Table II with time-invariant LEKs only. We assume the same 
set-up as in Section IV- A. We shall also assume thatp | 2n + l 
for the same reason as that in Theorem 4. 

Category 1 (Min-Cut between S2-D1 is equal to 0): This 
impUes that M21 = 0. Let the elements of Vi be given by 

[Vi],j^eij, j€{l,2,..,2n+l}, i€{l,2,..,n+l}, (28) 
where, 6ij is a variable that takes values from F^. Also, let 

V2 = M^s^MisViA and V3 = Ms^^MnViB, (29) 

where, the elements of the matrices A and B, of size 
(n + 1) X n, are given by [A]ij = Uij and [B]ij = bij 
respectively (aij and bij are variables that take values from 
¥q). The following theorem provides the conditions under 
which network alignment can be achieved. 

Theorem 6: For an acyclic three-source three-destination 
unicast network with delays, when the min-cut between S2-D1 
is equal to and the min-cut between the other sources and 
destinations are not zero, the input symbols X{"'~^^, Xl^ , and 
X3" can be exactly recovered at the sinks-1, 2, and 3 from 
the output symbols Y^"^^^, and 1^3^"+^ respectively, 

if 

Rank[Vi M^^ hh-yV-^ = Rank[Mi-^^M22V2 Vi] = 2n -M, 
Rank [M^^ Ah^Vs Vi]=2n + 1. 

Proof: Proof is given in Appendix J. ■ 
When the conditions of the above Theorem are satisfied, 
throughputs of ^til, and are achieved for 

the source-destination parrs S\ — Di, S2 — D2, and — -D3 
respectively. Hence, as n ^ 00, a throughput of 1/2 is 
achieved for every source-destination pair. 



Category 2 (Min-Cut between S2-D1, S3-D1 and S1-D2 
are equal to Oj: This implies that M21 = 0, M31 = and 
M12 = 0. Let the choice of V\ and V2 be the same as in (28) 
and (29) respectively, and choose the elements of V3 as 

[V3\ij=6ij, iG {l,2,..,2n+l}, iG{l,2,..,n}), (30) 

where 5ij is a variable that takes values from F,. The fol- 
lowing theorem provides the conditions under which network 
alignment can be achieved. 

Theorem 7: For an acyclic three-source three-destination 
unicast network with delays, when the min-cut between ^2- 
Di, S3-D1 and S1-D2 are equal to and the min-cut between 
the other sources and destinations are not zero, the input 
symbols X^", and X^" can be exactly recovered at 

the sinks-1, 2, and 3 from the output symbols Y^''+'^, Y^''+'^, 
and y3^"+^ respectively, if 

Rank[M32M22V2 V3] = 2n, Rank[M3^^Mi3Fi V3] =2n + l. 

Proof: Proof is given in Appendix K. ■ 

When the conditions of the above Theorem is satisfied, 
throughputs of and are achieved for 

the source-destination parrs Si — D\, S2 — D2, and ^3 — Ds 
respectively. Hence, as r? ^ 00, a throughput of 1/2 is 
achieved for every source-destination pair. 

Category 3 (Min-Cut between S3-D1, S1-D2 and S2-D3 
are equal to 0): This implies that M31 = 0, M12 = and 
M23 = 0. Let the choice of Vi be the same as in (28) and 
define the elements of V2 and V3 as 

[V2]ij=jij, [V3]ij=Sij, i € {l,2..,2n + € {l,2..,n} (31) 

where, 7^^ and Sij are variables that take values from F,. 

Theorem 8: For an acyclic three-source three-destination 
unicast network with delays, when the min-cut between ^3- 
Di, S1-D2 and S2-D3 are equal to and the min-cut between 
the other sources and destinations are not zero, the input 
symbols Xj""*"^, Xj", and X3" can be exactly recovered at 
the sinks-1, 2, and 3 from the output symbols Y^''^^, Y^'"+^, 
and ¥3^"+^ respectively, if 

Rank[Vi Mr/Mai V2] = RanklMg^MiaV^i V3] =2n-\-l, 
Rank[M3^^M22V2 V3] = 2n. 

Proof: Proof is given in Appendix L. ■ 
When the conditions of the above Theorem is satisfied, 
throughputs of (2n+i-) , and ^5;^ are achieved for 

the source-destination pairs Si — Di, S2 — D2, and 5*3 — D3 
respectively. Hence, as n — >■ 00, a throughput of 1/2 is 
achieved for every source-destination pair. 

Category 4 (Min-Cut between S3-D1, S3-D2, S1-D3 and 
S2-D3 are equal to 0): This imphes that M31 = 0, M32 = 
0, Mis = and M23 = 0. Here, we can achieve a sum- 
throughput of 2. Since, Ds is not facing any interference, we 
can take independent input symbols of source-3, i.e., Xg^""*"^, 



TABLE II 

Various possibilities of Min-Cuts between Source-? and DESTiNATioN-j {(i, j) e {l, 2, 3}|i ^ j) 





Min-Cut between Soiirce-i and Destination-j 


Category No. 


S2 - Di 


S3-D1 


Si - D2 


S3 - D2 


Si - D3 


S2 - D3 


1. 





> 1 


> 1 


> 1 


> 1 


> 1 


2. 











> 1 


> 1 


> 1 


3. 


> 1 








> 1 


> 1 





4. 


> 1 





> 1 












to be of size (2n+ 1) x (2n+ 1) and V3 to be an identity matrix Acknowledgement 

of size (2n+ 1) x (2n+ 1). The independent symbols, X'l"^^ 
and are column vectors of sizes (n + 1) x 1 and n x 1 
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Theorem 9: For an acyclic three-source three-destination g 5 Rajan 
unicast network with delays, when the min-cut between S^- 
Di, S3-D2, S1-D3, and S2-D3 are equal to and the min- 
cut between the other sources and destinations are not zero, 
the input symbols X'l"^^, X'^ , and X'^"^^^ can be exactly 

recovered at the sinks-1, 2 and 3 from the output symbols References 
y^2„+i^ y|"+i and Fg^^+i respectively, if 



Rank[Vi Mf/ MziVi] = Rank[Mr2^^22V2 Vi] = 2n -h 1. 

Proof: Proof is given in Appendix M. ■ 
When the conditions of the above Theorem is satisfied, 
throughputs of j^^y and g+1) are achieved for 

the source-destination pairs S\ — D\, S2 — D2, and 6*3 — -D3 
respectively. Hence, as n — > cxd, a throughput of 1/2 is 
achieved for Si — Di, S2 — D2 and a throughput of 1 is 
easily achieved for 5*3 — D?,. 

Remark 8: In all the above four categories, the choices of 
Vi, V2 and V3 were such that we could atleast ensure that Vx, 
V2 and V3 were full-rank, which were necessary to satisfy the 
network-alignment conditions. 

Remark 9: In Category 4, a sum-throughput of close to 2 
is achieved as n — > 00. For acyclic networks without delays, 
this can be easily achieved by time-sharing the network equally 
between source-1 and source-2. But it is not clear how such a 
sum-throughput can be achieved for arbitrary acyclic networks 
with delays whereas, our method provides a scheme that can 
achieve it. 

V. Discussion 

Though the transform method was originally claimed to be 
applicable for ayclic networks having M(D) whose elements 
are only polynomial functions in D, it can also be applied to 
networks having M (D) whose elements are rational functions 
in £) by multiplying by the LCM of all the denominators of 
the rational functions, at all the sinks. This gives a finite dmax- 
The same applies to cylic networks too. 

Network ahgnment for the three source-three destination 
unicast network with delays, discussed in this paper, can be 
extended to the case where each source-destination pair has 
a min-cut greater than one. We are currently working on it. 
An interesting dierction of future research is extending the 
network alignment to the case of arbitary number of sources 
and destinations with arbitrary message demands. 
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Appendix A 
Proof of Theorem 1 



Appendix C 
Proof of Lemma 2 



Proof: 



A 



The inverse of the matrix F is given by 

1 



\i=0 



(32) 



-1 =n-i 



1 1 

1 
1 



a 



1 a-2("-i) 



a 
a 



-in-l) 
-2(n-l) 



-(n-l)(n-l) 



Note that exists [9]. Now, can also be written as 
Qn = FiSiI^ (i.e. Kronecker product of F and 1^). Similarly 
Qu = F iSi lu- From (32), we have 

AQ^ = Qi/A.. 

Now, detiQ^) = [det{F)Y[det{I^)]-" ^ and Q-^ = (g) 
So, 

A = QAQ-\ 

Hence the theorem is proved. ■ 

Appendix B 
Proof of Theorem 2 

Proof: If both the conditions are satisfied after 
the assignment of values to e, then sink-j can invert 
[M^^Jihi) Ml'^]{h,) ■■■ M^^J,{k,,)] matrix and decode the 
required input symbols without any interference. 

If Condition 1) is not satisfied, then sink receives superposi- 
tion of required information and interference from other input 
symbols, which it cannot distinguish. 

If Condition 2) is not satisfied, then sink cannot invert 
the matrix [M\^]{h,) M^^Hr^ ■■■ M^']{K,)] which is 
necessary for decoding the input symbols. ■ 



Proof: Following the terminology developed so far, for 
some n » dmax and for each < t < n — 1, let 



X. 



it) 



Then, by (6), (14) and the structure of the M-^^ matrices, 
we have for < t < n — 1, 



, d=0 



where m|''^ is ai/jXfj, matrix over Fg (considered as a subfield 



of ¥qa such that 



(34) 



d=0 



We define a collection of ring homomorphisms 4>t '■ 
¥g{D) F,a for < i < n - 1, given by M^) = a*- 
For some matrix P{D) over Fq(_D), we also define cf)t{P{D)) 
to be equal to the matrix P with elements in ¥qa that are 
the ^j-images of the corresponding elements of P{D). Then, 
from (33) and (34), we have 



y.("-i-*) = </,,(M,(£>))X 



(n-l-t) 



(35) 



for < t < n — 1. Clearly, the zero-interference conditions 
satisfied in the Mj{D) matrices continue to hold in the 
(t)t{Mj{D)) matrices, for any Q <t <n—l and for any sink- 
j. Having satisfied the zero-interference conditions, to recover 
the source processes demanded by each sink-j at time instant 
n — l — t, the invertibility conditions also have to be satisfied, 
i.e.. 



l[det{MM'j{D)))^0, 



(36) 



where Mj{D) is the square submatrix of Mj{D) indicating 
the source processes that are demanded by sink-j. But then, 
we have 



det {(t>t{M'j{D))) = Mdet{M'.{D))) (37) 



and thus 



n det (MM'jiD))) = [] (det(Mj (£>))) 

l[det{M;iD)) 

= M.fm 



where f{D) is as defined in (15). Clearly, /(a*) ^ implies 
that (36) is satisfied and the source processes demanded at 
each sink can be recovered at time instant n — 1 — t in the 
transform approach. Similarly, if the sink demands are satisfied 
at time instant n — 1 — t in the transform approach, clearly 
we must have /(a*) 0. This holds for < t < n — 1, thus 
proving the lemma. ■ 

Appendix D 
Proof of Theorem 3 

Proof: If part: 

Let Fpm be the field over which the feasible network code 
has been obtained for {Q,C). Consider the polynomial f(D) 
(given by (15)) with coefficients from F^m. Let F^m' be 
the sphtting field of this polynomial, i.e., a suitable smallest 
extension field of Fpm in which f{D) spUts into linear factors. 

Let 



b=k 



P" 



1 



5=1 



where each pb is some prime and rrift is some positive integer. 

By Lemma 2, the choice of a to be used for the DFT oper- 
ations should be such that /(a*) ^ 0, for any < t < n — 1. 
We now show that such an a exists and can be chosen. 
Let F ,„// be an extension field of F^™/ . Clearly, 

p'" ^ 1 j I (p"^ ~ ^) ■ ^o^^^^'"' further demand that 
is such that 



6=fe 



lipfUPc 

6=1 c=l 



(38) 



where each pc is some prime and m'^ and m" are some positive 
integer such that pb ^ Pc ior 1 < b < k and 1 < c < k'. Note 
that m'^ > nib for 1 < b < k. Such extensions of Fp„," can 
indeed be obtained. For example, F^n," can be considered to 
be the smallest field which contains F ™' and F^™ , m being 
some positive integer coprime with m . Then clearly F^m" is 
such that (38) holds. 

Following the notations of Section III, we now pick a G 
Fpm" (where m" satisfies (38)) such that the following con- 
dition holds 

• The cyclic subgroup {l, a, a"~^} of F^^" \ {0} with 
order n{n > 1) is such that n and nb^iP™'' ^® coprime. 
Such an a can be obtained by choosing a from the subgroup 
of Fp,„" \ {0} with n = Ylt^i P^" elements. We now claim 
that using such an a for the DFT will result in a feasible 
transform network code for {Q,C). The proof is as follows. 

We first note that the zero-interference conditions are satis- 
fied irrespective of the choice of a in the DFT operations. As 
for the invertibiUty conditions, by Lemma 2, it is clear that 
as long as /(a*) 7^ for 0<t<n — 1, we have a feasible 
transform network code for {G,C). Suppose /(a*) = for 
some 1 < t < n — 1. Let be the order of a*, i.e. the number 
of elements in the cyclic group generated by a*. Then ni\n and 

also 'it|n6=iP(r'' ^ ^p"'' ^ °f /(-^)- However 



this leads to a contradiction as n shares no common prime 
factor with ri6=i-Pr' • Thus no a*, I < t < n — 1, can be a 
zero of f{D). This, coupled with the given fact that /(I) 7^ 0, 
proves the claim and hence the if part of the theorem. 

Only If part: 

Let Fq be the field over which a feasible transform network 
code has been defined for {Q,C), i.e., there exists a choice of 
LEKs and for a from Ff, using which the zero-interference 
and the invertibility constraints have been satisfied in the 
transform domain. Note that a choice for the LEKs impUes 
that the matrices Mj(D) given by (6) are well defined. We 
will now prove that the invertibility and the zero-interference 
constraints also hold in these Mj{D) matrices for all sinks, 
i.e. for 1 < J < r. 

We first prove the invertibihty conditions. Towards that end, 
let M^"~^' be defined as the uj x p, transfer matrix at time 
instant n — 1 from all the sources to sink-j in the transform 
approach, i.e.. 



(n-l) 



■1) 



...M. 



(n-l) 



(39) 



By the structure of the M^^ matrices, we have Afj" = 
^^jd) ^ Mj{D)\d=i. Let Af^-*""^' be the submatrix 
of Afj"^^^ which is known to be invertible, as it is given that 
the invertibility conditions for the transform network code are 
all satisfied. 

The invertibility conditions for sink-j of the usual (non- 
transform) network code for {Q,C) demand a suitable sub- 
matrix M'j{D) of the matrix Mj{D) to be invertible. Note 



however that M' {D)\d= 



M. 



(n-l) 



by (39). Therefore, we 



have det (^m/" = det {M'j{D)\d=i) + 0. As in (37), we 

have det {M]{D)) = det {M'^{D)\d=i) 7^ 0. Therefore, 
det {M^{D)) ^ 0, i.e., det {M^{D)) is a non-zero polynomial 
in D. Because the choice of the sink was arbitrary, it is clear 
that the invertibility conditions hold for each sink in the usual 
network code for {G, C). By (15), we also have (D-l) | f{D). 

We now prove the zero-interference conditions. The zero- 
interference conditions in the transform domain can be inter- 
preted as follows. Having ordered the input processes at the 
source-i, suppose the sink-j does not demand the A;*^ process 
from the source-i. Then the matrix Mij is such that fc*'* column 
of M-j^ is an all-zero column for all < t < n — 1. To prove 
that the zero-interference conditions continue to hold in the 
usual network code for {G,C), we must then prove that for 
each source-i, each particular sink-j and each k (such that the 
fc*'' input process at source-i is not demanded at sink-j, the 
k*^ columns of M-f^ matrices are all-zero for < c? < dmax 



where M^f\0 < d < dr, 



are matrices such that 



d=0 



This is seen by observing the structure of the Mij matrix, 
which is defined by (1 1). Using Theorem 1 and with /?„ = a°, 
we have (40) (shown at the top of the next page). Comparing 



Mij = Q^.MijQ^l 



2j 



an-l-tjQ-(*) 



on-l T 



En— 1 yOTI 
t=0 Pi 

on-l-t iw-(t) on-l-t o-l iw- 



(t) 



M. 





(n-2) 



••• 

••• 







M. 



(0) 



En-l o-{n-l) ii^fCt) 
t=0 Pn-l-t ^^^ii 

En-l on-l-t o-(n-l) ;,-V(t) 
t=0 "l Pn-l-t "'^ij 



1 )\>(«) 



(„_i) 

Pn-l 



(40) 



the submatrices of Mij from (11) and (40), we see that if 
the fc*'' column of the M-'*'' matrices is all-zero for all < 

t < n — 1, then the fc*'' columns of M'>f^ matrices are all- 
zero for < d < d-max- As the choice of source-i and sink-j 
are arbitrary, it is clear that the zero-interference conditions 
continue to hold in the Mij (D) matrices for all 1 < i < s and 
^ < j < r. This proves the only if part of the theorem and 
hence the theorem is proved. ■ 

Appendix E 
Proof of Lemma 3 

Proof: Consider Mij as defined in (11) which is a 
circulant matrix of size (2n + 1) x (2n + 1). Note that the 
diagonal elements of Mij, i.e., M-j^ {k e 0, 1, .., 2n), are the 
eigen values of the matrix M^j. Also, note that the eigen values 
are equal to (2n + l)-point finite-field DFT of the first row 
of Mij. Since, the min-cut from source-i to sink-j is equal to 
1, by Menger's Theorem, there exists exactly one link-disjoint 
directed path from source- z to sink-j. Let such a directed path 
consist of links ei, 62, .., ef. Now, we can assign the values 

Q!l,ei = 1. /3ei,ei+i = 1 (« € {1,2, 1}) , Ee^,! = 1 

and assign values of to all the other LEKs. By, such an 
assignment of values to the LEKs, exactly one among M^^, 

Mj^p, .., M^j""'"'^ is equal to 1. This implies that all the eigen 
values of Mij are non-zero. Hence, the diagonal elements of 
Mij are non-zero polynomials in e and so is its determinant. 

■ 

Appendix F 
Proof of Lemma 4 

Proof: If part: Euler's theorem [14] states that if two 
positive integers a and b are coprime then, b divides a'^^''^ — 1 
where (p represents the Euler's totient function. If 2n -|- 1 <p 
then, 2n + 1 and p are coprime. If 2n + 1 > p then, p and 
2n+ 1 are coprime iff p does not divide 2n+ 1. Hence, by 
Euler's theorem, 2n + l|p'^(2"+i) - 1 if p f 2n -M. Thus if 
p\2n+l then, 2n+ 1 |p™ - 1, for all m such that 4>{'2n+l)\m. 
Only If part: If 2n+ 1 divides — 1 for some positive integer 
m then, — 1 = r(2n -I- 1) for some positive integer r. So, 
p-m _ ^2fi ^ 1)7. = 1 which means that p and 2n -|- 1 must be 
coprime. Since p is prime, p\2n+l. ■ 



Appendix G 
Proof of Theorem 4 

Proof: To exactly recover X^" and X^" at the 

sinks-1, 2 and 3 respectively, it is sufficient that the following 
network aUgnment conditions are satisfied. 

M21V2 - M31V3 (41) 

Mz^Vz C Mi2Vi (42) 

M23.V2 C MizVi (43) 
Rank[Miiyi M21V2] = 2n+l 

^ Rank[l/i Mf/ ^^21^2] = 2n+l (44) 
Rank[M22V2 M12V1] = 2n -h 1 

^ Rank[Mi"2^M22V'2 Vi] =2n+l (45) 
RankfMssl/s M13V1] =2n+l 

^ Rank[M{^^M33F3 Vi] =2n+l (46) 

Note that from Lemma 3, inverse of Mij V G {1,2,3} 
is well-defined. It is easily seen that the choice of Vi, V2 ,and 
V3 in (19)-(21) satisfy the conditions (41)-(43). Suppose that 
(44)-(46) are satisfied. Let 

fi{e) = det{[Vi M-'M2iV2]) 
f2{s) = det{[Mr2M22V2 Vi]) 
/3(e) = det([M^,'M33V3 Vi]) 
h{e) = n det{Mij) 

(iJ)G{l,2,3} 

/(£)=n/ife). 

i=l 

Since /i(e), /2(e) and /3(e) are non-zero polynomials in e, 
/(e) is also a non-zero polynomial in e. Hence, by Lemma 
1 in [3], for a sufficiently large field size, there exists an 
assignment of values to s such that the network alignment 
conditions are satisfied. Since p | 2n + 1, by Lemma 4, for a 
sufficiently large m (in particular, m such that (j){2n + l)\m 
where ^ represents the Euler's totient function), there exists 
an assignment of values to e such that the network alignment 
conditions are satisfied. Hence, the theorem is proved. ■ 



Appendix H 
Proof of Lemma 5 

Proof: If we assign £(^'''"f^) = gC-^max+i) — _ = 
£(2") = £^ jvf^_^. in (25) becomes a circulant matrix and hence 
can be diagonalized as shown in Theorem 1 . Further, in lemma 
3 we proved that the determinant of the diagonahzed matrix 
is a non-zero polynomial in e. So, the determinant of the 
circulant matrix is also a non-zero polynomial in e. Hence, 
the determinant of Mij is a non-zero polynomial in e'. ■ 

Appendix I 
Proof of Theorem 5 

Proof: To exactly recover X["''^^, Xl^ and X'^ at the 
sinks-1, 2 and 3 respectively, it is sufficient that the following 
network alignment conditions are satisfied. 

Span(M2i ^2 ) = Span(M3i ^3 ) (47) 

Span(M32V3) C Span(Mi2Vi) (48) 

Span(A/r23V'2) C Span(Mi3Fi) (49) 
Rank[MiiFi M21F2] = 2n + 1 

Rank[t/i M1VM21F2] = 2n + 1 (50) 
Rank[M22V"2 Afi2Vi] = 2n + 1 

<!4> Rank[Afi"2^Af22V"2 T4]=2n-M (51) 
Rank[M33y3 Mi3yi] = 2n -h 1 

<S4> Rank[M{^iM33y3 Fi] = 2n -M (52) 

The choice of Vi, V2 and V3 ensures that the conditions (48) 
and (49) are satisfied. To satisfy (47), we have to ensure that 

M^^M-ixVi = Ms2^Mi2ViBC 
^ M.^^H{2iM^/Mi3ViA = M^^Mi2ViBC 
^ TiViA = ViBC (53) 

is satisfied. In order to satisfy (53), every element of TiViA 
must be equal to every element of ViBC, i.e., 

Qij = 0, V i e {1, 2, .., 2n + 1}, j G {1, 2, .., n}. 

Hence, to satisfy (50)-(53) we need to find an assignment to 
the variables - 0, e', a and b, such that / ^ and gij = 0, 
V i e {1,2, ..,2n-|-l}, j e {1,2, ..,n}. This means that there 
must exist an assignment such that ^ and g^^^^ = 0. 

After the assignment to the variables, we require that 
and i?!^'^^ ^ as dividing by zero is prohibited. In order to for- 
mulate this as an algebraic problem, introduce a new variable 
S and consider the polynomial ^1 — (J/f"*") J]^^^ ^.^ g^j^^^ . 
From Weak Nullstellensatz [15], an assignment to the variables 
- ^) e') o,, b and c exist such that p]'"'^^ = 0, for all 
and (1 - n(i,^) g-f^) = iff 1 does not belong 

to the ideal generated by the polynomials g^^^^ for all 
and(l-5/(-)/Wn(,,)fff^). 



Appendix J 
Proof of Theorem 6 

Proof: Let 6 = {Oij V {i, j)}, a = {atj V and b = 

{bij V («,.?)}. To exactly recover X'^"^^, X'^ and X3" at the 
sinks-1, 2 and 3 respectively, it is sufficient that the following 
network alignment conditions are satisfied. 

S'pan(M32V^3) C Span{Mi2Vi) (54) 
Span{M23V2) C Span{MxiVi) (55) 
Rank[Miiyi M31V3] =2n+l 
^ Rank[Fi M^^M^iVs] = 2n+l (56) 

Rank[M22V2 M12V1] = 2n + 1 
Rank[Mi2^Af22V2 yi] = 2n+l (57) 
Rank[M33F3 M13V1] = 2n + 1 
4^ Rank[A^{3iAf33V:3 Vi]=2n + l (58) 

It is easily seen that the choice of V2 and V3, in (29), satisfy 
the conditions (54) and (55). Suppose, (56)-(58) are satisfied. 
Now, let 

fiis,0,a)=deti[M-'M22V2 V^]) 
f2{e,e,b)^det{[Vi M-.'li'hiVa]) 

/4(£) = n det{Mij) 

(i,i)e{l,2,3}|(«j)5^(2,l) 

Let f{e,e,a),b) = f i{e, 0, a) Uk)fik^ IS- Since, 
/i (£7 ci), ./2 (e, 0, b) and fs (e, 6, b) are non-zero polynomials, 
f{£,0,Q:,b) is also a non-zero polynomial in e. Hence, by 
Lemma 1 in [3] and Lemma 4, for a sufficiently large field 
size, there exists an assignment of values to variables e, such 
that the network ahgnment conditions are satisfied. Hence, the 
theorem is proved. ■ 

Appendix K 
Proof of Theorem 7 

Proof: Let S ~ {Sij V To exactly recover X["'^^, 

Xi^ and X3" at the sinks-1, 2 and 3 respectively, it is 
sufficient that the following network alignment conditions are 
satisfied. 

Span(A^23V2) C Span(Afi3l4) (59) 

Rank[Afiiyi] =n-hl (60) 

Rank[A^22V"2 M32V3] = 2n 
^ Rank[Af3"^^M22V'2 ^3] = 2n (61) 

Rank[Afi3 Fi Mz^Vi] = 2n+l 
<^ Rank[if3-3^ifi3l/i ^3] = 2n+l (62) 

It is easily seen that the choice of V2 as in (29), satisfies the 
condition in (59). Since Vi is full-rank and Afn is invertible 



(from Lemma 3), (60) is also satisfied. Suppose (61) and 
(62) are satisfied. Let fi{s,9) denote the determinants of 
all the n X n sub-matrices of MnVi. Also, let f2{§.,d,5,a) 
denote the determinants of all the 2n x 2n sub-matrices of 
[M^2^M22V2 Vs]. Now, define 

U{e) = n det{Mij). 

(»j)G{l,2,3}|(iJ)#{(2,l),(3,l),(l,2)} 

Let f{s,e,6,a) = h{e,e)h{e,e,5,a)h{e,e_A)fi{e)- Since, 
/i(£)^)> f2{§.,0,d,a} and fz{£,0_,5,a) are non-zero polyno- 
mials, f{e,d_,5_,a) is also a non-zero polynomial. Hence, by 
Lemma 1 in [3] and Lemma 4, for a sufficiently large field 
size, there exists an assignment of values to variables e, 6, 5 
and a such that the network alignment conditions are satisfied. 
Hence, the theorem is proved. ■ 

Appendix L 
Proof of Theorem 8 



X^" and X^" at the 



Proof: To exactly recover X[ 
sinks-1, 2 and 3 respectively, it is sufficient that the following 
network alignment conditions are satisfied. 





Rank[MiiVi M21V2] 


= 2n + 1 






Rmk[Vi Mfi^MaiFa] 


= 2n+l 


(63) 




Rank[A^22^2 M32V3] 


= 2n 






Rmk[M^2^M22V2 V3] 


= 2n 


(64) 




RankiMisFi MsaVa] 


= 2n+l 




<«4> 


Raiik[M33^Mi3Fi V3] 


= 2n+l 


(65) 



Suppose that (63)-(65) are satisfied. With 7 = {-jij V 

let /i (e, 7, S) denote the product of determinants of all the 

2n X 2n sub-matrices of [M^2^M22V2 V3]. Also, let 

/2(£,S,7) = det{[Vi M1VM21F2]) 
f3ie,e,S)=det{[M^,'M,3Vi V3]) 
/4(£) = n det{Mij). 

(i,j)e{l,2,3}|(i,j)#{(3,l),(l,2),(2,3)} 

Let f{e,e,2,S) = h{s,2,S)f2ie,e,2)f3{e,0,6)U{e). 
Since, /i (e, 7, 5 ), /2 (e, 0, 7) and /s {s, 0, S) are non-zero poly- 
nomials, f{e, 9, 7, S) is also a non-zero polynomial. Hence, by 
Lemma 1 in [3] and Lemma 4, for a sufficiently large field 
size, there exists an assignment of values to the variables s, 
6, 7 and 6 such that the network aUgnment conditions are 
satisfied. Hence, the theorem is proved. ■ 

Appendix M 
Proof of Theorem 9 



network aUgrmient conditions are satisfied. 

RankfMas] = 2n -h 1 (66) 
RankfMiiVi M2iV^2] =2n+l 

^ Rmk[Vi Mf/MaiFa] = 2n-h 1 (67) 
Rank[M22V2 ^12^1] = 2n + 1 

4^ Rank[Mf/M22F2 yi]=2n+l. (68) 

Since M33 is invertible (from Lemma 3), (66) is satisfied. 
Suppose that (67) and (68) are satisfied. Let 

fi{s) = det{[M3a]) 

f2{e,0n)=det{[Vi Mu'M2iV2\) 

h{eA,j) = det{[M^2^M22V2 Fi]) 

h{e) = n det{Mij). 



With 



7 



(i,j)e{l,2,3}|(i,j)5^{(3,l),(3,2),(l,3),(2,3)} 

= V (i,.?)}, let /(e,£,7) 



/i(£)/2(£,^,7)/3(£,^,7)/4(£). Since, /i(e), /2(£,^,7) 
and f3{s,6,j) are non-zero polynomials, f{s,0,j) is also a 
non-zero polynomial. Hence, by Lemma 1 in [3] and Lemma 
4, for a sufficiently large field size, there exists an assignment 
of values to the variables s, 6 and 7 such that the network 
alignment conditions are satisfied. Hence, the theorem is 
proved. ■ 



Proof: To exactly recover X'^"^^, X2" and X^^"^^ at the 
sinks-1, 2 and 3 respectively, it is sufficient that the following 



